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Abstract 


In  an  experiment  which  was  conducted  to  determine  the 
response  of  an  isolated  rotor  to  an  inlet  distortion, 
a triaxial  hot-film  probe  was  used  to  make  velocity  and 
flow  angle  measurements  near  the  exit  plane  of  the  rotor. 
The  velocity  sensors  were  three  mutually  orthogonal  film- 
type  sensors  operated  in  the  constant  temperature  mode. 

This  probe  has  been  calibrated  and  procedures  have  been 
developed  to  extract  velocity  components  from  the  anemo- 
meter output  voltages.  It  has  been  determined  that  the 
proximity  of  the  probe  shaft  results  in  severe  distortion 
of  the  flow  over  the  velocity  sensors.  The  process  which 
has  been  developed  for  computing  velocities  yields  satis- 
factory results  in  test  cases  and  is  expected  to  be  suit- 
able for  reducing  the  data  from  limited  number  of  flow 
conditions  for  which  the  probe  was  used  in  the  present 
experiment.  However,  it  is  considered  to  be  too  involved 
and  too  susceptible  to  error  to  be  used  in  experiments 
where  there  is  either  a vast  amount  of  data  to  be  processed 
or  there  is  no  means  of  spot -checking  the  results.  It  is 

therefore  considered,  mandat-cry-  Cha-t-wn  •improved''  probe 

geometry  be  devised  which  does  not  exhibit  the  degree  of 
interference  observed  with  the  probe  in  question  before 
undertaking  to  make  extensive  wake  measurements. 
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Introduction 

An  inlet  distortion  experiment  has  been  conducted  in 
United  Technologies  Research  Center's  (UTRC)  Large  Scale 
Rotating  Rig  (LSRR) . (Ref.  1)  During  the  course  of  this 
experiment  some  runs  were  conducted  with  a triaxial  hot- 
film  probe  mounted  in  the  rotating  reference  frame  near 
the  rotor  exit  plane.  This  probe  was  used  to  make  velocity 
and  flow  angle  measurements  in  the  unsteady  flow  field  of 
the  rotor  as  it  passed  through  a distorted  inlet  flow.  The 
three  mutually  orthogonal  velocity  sensors  mounted  on  the 
probe  were  operated  in  the  constant  temperature  mode  by 
anemometer  circuits  developed  in-house  especially  for  this 
experiment . 

Before  the  data  obtained  from  such  a probe  may  be  used 
to  determine  velocity  and  flow  angles,  it  must  be  calibrated 
to  determine  the  response  of  each  sensor  to  normal  velocity 
and  to  ascertain  the  manner  in  which  they  respond  to  trans- 
verse flows.  The  measurements  necessary  to  accomplish  this 
calibration  were  carried  out  in  North  Carolina  State 
University's  (NCSU)  low  speed  wind  tunnel. 

Equipment 

The  probe  used  in  the  LSRR  inlet  distortion  experiment 
is  a Thermo -Systems , Inc.  1296M  triaxial  hot-film  probe. 

A total  pressure  transducer  is  attached  to  the  support 
shaft  resulting  in  the  configuration  shown  in  Fig.  1.  Two 
axes  systems  are  used  in  calibrating  the  triaxial  hot-film 
probe.  One  system  is  the  axial -radial -c ircumferential  axis 
system  of  the  rotor  while  the  other  is  chosen  with  the  three 
axes  parallel  to  the  three  mutually  orthogonal  velocity 
sensors  which  make  up  the  probe.  These  two  systems  are 
shown  in  Fig.  2 omitting  the  total  pressure  probe  for 
clarity.  The  velocity  components  in  the  two  systems  may 
be  related  by  two  sets  of  transformation  equations.  The 
first  set  gives  velocities  in  the  rotor  reference  frame, 
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v^,  in  terms  of  the  sensor  frame  velocities,  u^. 


Vi  = ai i u i 

+ 

a2 iu2  + 

a3  iu3 

(1) 

v2  = ai 2ui 

+ 

a 2 2 u 2 + 

a 3 2U  3 

(2) 

v3  + ajjUi 

+ 

a2  3u2  + 

a3  3U3 

(3) 

The  second  set 

of 

equations  is  the 

reverse  transforma 

ui  = ai 1V1 

+ 

ai2v2  + 

a 1 3 V3 

(4) 

u2  + a2 i Vi 

+ 

a22v2  + 

a2  3V3 

(5) 

u 3 = a 3 i Vi 

+ 

a3  2v2  + 

a 3 3 V 3 

(6) 

The  coefficients,  » used  in  these  equations  are  as  follows 


aii  - VTTS 
a 1 2 = 0 


a2 i = -/1/6 
a2  2 = -/T7T 


a3 i = -/1/6 

a32  = /T7I 


a i 3 = /T/3 


a2  3 = /T73 


a3  3 = /173 


Each  of  the  sensors  is  operated  in  the  constant  resis- 
tance (constant  temperature)  mode  by  a circuit  developed 
in-house  especially  for  this  experiment.  A commercially 
produced  unit  was  first  tried  but  it  failed  under  the  centri- 
fugal loads  encountered  during  operation  of  the  LSRR.  The 
present  circuit  is  compact  and  quite  durable  when  encapsu- 
lated in  epoxy.  The  simplicity  of  the  ciruit  is  evident 
in  the  schematic  diagram  shown  in  Figure  3.  This  circuit 
has  been  subsequently  modified  to  make  it  more  versatile  and 
improve  performance  slightly  but  the  original  circuit,  as 
shown,  proved  to  be  quite  adequate  for  the  present  application. 
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The  probe  was  calibrated  in  the  NCSU  Low  Speed  Wind 
Tunnel.  This  tunnel  is  of  conventional  closed  loop  design 
with  the  fan  located  downstream  of  the  test  section  which 
is  32  inches  high  and  45  inches  wide.  Turning  vanes  are 
used  in  the  corners  and  flow  straightening  screens  are  lo- 
cated up-stream  of  the  test  section  to  provide  a uniform  flow. 
A probe  support  fixture  was  designed  which  allows  the  probe 
to  be  pitched  and  yawed  while  maintaining  the  sensors  in 
precisely  the  same  location  thus  allowing  any  flow  non- 
uniformities which  might  still  exist  to  be  neglected.  This 
fixture,  which  is  shown  in  Fig.  4,  may  be  pitched  + 20°  in 
5°  increments  and  yawed  continuously  through  a full  360° 
although  angle  readout  is  available  for  only  + 90°  The 
sign  conventions  used  for  the  pitch  and  yaw  angles,  0 and  <j>, 
respectively,  are  shown  in  Fig.  5.  These  angles  and  the 
total  velocity  vector,  V,  may  be  used  to  calculate  the  velo- 
cities in  the  rotor  reference  frame. 


Vi 

= V 

sin 

0 

(7) 

v2 

= V 

cos 

0 

sin  <J) 

(8) 

V3 

= V 

cos 

0 

COS  <J) 

(9) 

Procedure 

The  velocity  in  the  Low  Speed  Wind  Tunnel  is  measured 
by  a slant  manometer  connected  to  a pressure  port  in  the 
settling  chamber.  The  static  pressure  in  the  test  section 
is  atmospheric  as  the  test  section  is  vented  to  the  atmosphere. 
In  order  to  obtain  a more  accurate  velocity  measurement,  a 
pneumatic  total  pressure  probe  was  mounted  in  the  probe 
support  fixture  so  that  its  head  was  in  the  exact  position 
that  the  sensors  of  the  triaxial  probe  would  occupy  when  it 
was  inserted.  A static  pressure  port  in  the  test  section 
floor  was  also  monitored  and  the  tunnel  was  operated  over 
the  range  of  velocities  to  be  used  in  calibration  to  obtain 
a precise  calibration  of  the  velocity  at  the  probe  head  as 
a function  of  the  manometer  reading. 
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The  yawing  mechanism  of  the  probe  support  fixture  was 
adjusted  to  read  zero  with  the  total  pressure  probe  aligned 
directly  into  the  flow.  The  pitching  mechanism  needed  no 
adjustment  but  this  could  have  been  accomplished  by  shim- 
ming the  base  if  necessary. 

The  total  pressure  probe  was  removed  and  the  triaxial 
hot-film  probe  was  installed  into  the  calibration  fixture 
using  a flag  on  the  probe  shaft  to  align  it  with  the 
fixture.  The  probe  geometry  was  used  to  calculate  pitch 
and  yaw  angles  which  would  place  each  sensor  in  turn  per- 
pendicular to  the  flow.  Based  on  previous  experience,  it 
was  suspected  that  the  probe  shaft  was  deflecting  the  flow 
so  that  the  flow  was  not  actually  normal  to  the  sensor 
in  question  at  the  calculated  pitch  and  yaw  angle.  Con- 
sequently, with  the  tunnel  running,  the  probe  was  yawed 
slightly  while  monitoring  the  anemometer  output  voltage. 

The  flow  was  taken  to  be  normal  to  the  sensor  at  maximum 
output  voltage.  The  yaw  angle  at  peak  output  did  vary 
slightly  from  the  calculated  values  confirming  that  the 
probe  shaft  was  indeed  interfering  with  the  airflow  over 
the  sensors.  With  each  sensor  in  turn  placed  normal  to 
the  flow,  the  tunnel  was  operated  over  a wide  range  of 
velocities  and  the  various  flow  properties  and  the  output 
voltage  of  the  appropriate  anemometer  were  recorded.  These 
data  are  presented  in  Table  I and  are  plotted  in  Figs.  6 
through  8 for  sensors  1 through  3,  respectively.  It  was 
desired  to  fit  the  data  to  the  form  of  Eq . 10 

E2  - Eo  = K/pV  (Ts  - TJ  , (10) 

where  E is  the  anemometer  output  voltage,  p is  the  density, 
and  Tg  and  Tro  are  the  sensor  and  free  stream  temperatures, 
respectively.  The  constants  Ii  and  K arc  determined  from  the 
calibration  data.  No  single  value  of  E and  K was  found  to 
be  valid  over  the  entire  velocity  range;  consequently,  a 
least  squares  procedure  was  used  to  obtain  two  curve  fits 
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of  the  form  of  Eq . 10,  one  of  which  was  for  the  higher 
velocities  and  the  other,  for  the  lower  velocities.  The 
anemometer  voltage  at  the  intersection  of  these  two  curves 
was  calculated  and,  along  with  the  values  of  Eq  and  K,  is 
presented  in  Table  II  for  the  three  sensors.  During  the 
remainder  of  the  calibration,  if  the  output  voltage  of  the 
anemometer  was  greater  than  the  crossover  voltage  for  that 
channel,  the  upper  curve  fit  was  used.  On  the  other  hand, 
if  it  was  lower,  the  lower  curve  fit  was  used. 

After  determining  the  response  of  the  probe  sensors  to 
normal  flow,  it  was  necessary  to  investigate  the  behavior 
of  the  sensors  in  non-perpendicular  flows.  One  expression 
which  is  commonly  used  to  describe  the  response  of  a sensor 
to  transverse  flow  is: 

qeff  = un  + K*ut  ’ C11) 

where  qg££  is  the  effective  velocity  measured  by  the  sensor, 
u and  ut  are  the  normal  and  transverse  velocity  components, 
respectively,  and  K is  a constant  which  is  empirically  de- 
termined. Unfortunately,  this  expression  does  not  lend 
itself  to  situations  where  the  flow  may  be  blocked  by  the 
probe  shaft.  The  formulation  of  Eq.  12  was  developed  to 
include  these  blockage  effects  as  well  as  sensitivity  to 
transverse  flow. 

qeff  = un  ft0,  ^ , t12) 

where  qe££  and  u^  are  as  before  and  f(0,  <}>)  is  a weighting 
factor  of  order  unity  which  is  a function  of  the  pitch 
and  yaw  angles.  The  wind  tunnel  was  operated  at  a constant 
velocity  of  116  ft/sec  while  the  probe  was  pitched  and  yawed 
and  the  data  obtained  used  to  calculate  the  value  of  this 
function  for  each  sensor  at  each  combination  of  pitch  and 
yaw  angles.  The  values  thus  obtained  are  listed  in  Table  III. 


For  sensors  2 and  3,  the  values  are  generally  greater  than 
unity  with  the  low  values  occurring  only  where  the  probe 
shaft  might  be  expected  to  block  the  flow  (i.e.  the  flow 
is  nearly  normal  to  the  sensor)  and  the  higher  values 
occurring  where  there  is  significant  traverse  flow.  For 
sensor  1,  evaluation  of  the  function  yielded  values  much 
less  than  unity  over  the  entire  range  of  pitch  and  yaw 
angles.  This  behavior  indicates  that  something  is  blocking 
the  flow  over  sensor  1 at  all  times. 

The  values  of  the  function  at  each  value  of  0 were 
least  square  fitted  to  a fifth  order  polynomial  in  4> . 

The  coefficients  of  each  term  in  this  polynomial  were  then 
least  squares  fitted  to  a fifth  order  polynomial  in  0 
resulting  in  a 5 x 5 matrix  of  coefficients  for  each  sensor 
from  which  the  value  of  the  function,  f (0,  4>)  , at  any  inter- 
mediate value  of  0 and  4>  may  be  computed.  The  computer  codes 
used  to  reduce  the  tabulated  function  to  the  matrix  of 
coefficients  and  to  subsequently  reconstruct  the  function 
at  arbitrary  0 and  <f>  are  given  in  Appendix  A. 

In  order  to  use  the  calibration  data  to  obtain  the 
velocity  and  flow  angles  from  the  three  anemometer  output 
voltages,  one  first  must  calculate  an  effective  normal 
velocity  for  each  wire  by 

q = [ ^ ]2  (13) 

k/p"  (Ts  - TJ 

where  quantities  on  the  right  hand  side  of  the  equation 
are  the  same  as  for  Eq.  10.  Equation  12  may  be  written  for 
each  sensor  as  follows. 


P-2  = u + u 
t 1 2 


^ - U 

l 2 1 

p-2  = U 

I 3 1 


2 + u 2 


2 + U 2 
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In  these  equations,  f.  is  the  value  of  the  function  f(0,  <t>) 
for  the  ith  sensor  which  is  evaluated  at  the  value  of  0 and  <{> 
from  the  previous  iteration.  These  equations  may  be  solved 
quite  easily  for  the  squared  velocity  components,  u^2 . Since 
there  is  no  way  to  know  whether  u^  is  positive  or  negative 
from  these  equations,  the  proper  sign  must  be  known  from 
some  other  source.  In  the  present  experiment,  the  orien- 
tation of  the  probe  was  such  that  the  velocity  components 
may  always  be  assumed  to  be  positive.  Using  the  transfor- 
mation equations  1 through  3,  the  velocity  components  in  the 
rotor  reference  frame,  v^,  are  computed.  The  total  velo- 
city, V,  and  the  resultant  flow  angles,  and  <{>  , are  cal- 
culated by  Eqs . 17  through  19. 


V = /vi 2 + v2  2 + v3  2 

0r  = sin  (vi/V) 

_ i 

<J)  = sin  (v2/V  cos  0) 


(17) 

(18) 
(19) 


In  all  cases  tested,  the  method  overcorrects.  That  is,  if 
the  value  of  used  to  evaluate  the  functions  f(0,  <j>)  is 
correct  but  0 is  lower  than  the  correct  value,  0 will  be 
greater  than  the  correct  value  of  0.  If  the  value  of  0 is 
not  greatly  in  error,  <|>  will  be  essentially  the  same  as  • - 
the  assumed  value  of  4> . For  this  reason,  the  iterations 
were  uncoupled  and  under-relaxed.  Thus,  with  6 being  held 
constant,  <J>  is  iterated  upon  using 

4>n  + 1 = 0 . 8<f>n  + 0.2<J>rn  (20) 

where  the  superscripts  denote  the  iteration  level.  When  <J>n 
and  <J>  n have  converged  within  the  desired  tolerence  (normally 
2.5°),  the  value  of  <j>n+  is  set  by  Eq . 20  and  e is  updated 
by  Eq.  21. 
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en+1  = 0.8  en  + 0.2  e n * (21) 

r 

When  0 has  also  converged,  the  process  is  terminated. 

Note  that  when  the  iteration  has  been  terminated,  the 

correct  value  of  0 is  between  0n  and  0 n which  are  in 

r 

error  by  the  tolerence,  e.  Thus  choosing  the  computed 
value  of  0 to  be  the  average  of  0n  and  0rn  guarantees 
accuracy  within  e/2.  However,  due  to  the  overshoot  which 
appears  to  be  inherent  in  the  method,  allowing  the  computed 
value  of  0 to  be  defined  by  equation  21  gives  a lower 
average  error.  (The  maximum  possible  error  for  any  single 
calculation  is  0.8e  when  using  this  method  but  the  average 
error  appears  to  be  on  the  order  of  0.2e.)  A similar 
argument  may  be  made  for  the  evaluation  of  cp . 

As  a test  case,  the  anemometer  voltages  obtained  from 
a run  in  which  the  rotor  blades  were  removed  and  the  rotor 
was  run  very  slowly  in  order  to  survey  the  flow  field, 
were  inserted  into  Eq . 13  and  the  process  previously  out- 
lined was  followed.  Here,  all  three  velocity  components, 

Vi , v 2 , and  v3,  were  known.  There  should  be  no  radial  flow, 

Vi,  the  circumferential  velocity,  v2  , should  be  entirely  due 
to  the  rotation  of  the  rotor,  and  the  axial  velocity,  v3, 
should  match  that  measured  by  the  pneumatic  instrumentation 
in  the  stationary  reference  frame.  (The  pneumatic  instru- 
mentation used  to  measure  the  axial  velocity  was  normallv 

l • ••  «%»  • - • * • m • - • • • • • *.  ar  » ' « • 

used  only  to  set  the  flow  condition  and  is  not  very  accurate. 

It  normally  gave  an  axi->  velocity  somewhat  less  than  the 
true  value.  However,  it  was  ill  that  was  available  in  this 
one  case.)  All  three  calculated  velocity  components  were 
grossly  in  error. 

Since  there  were  previous  indications  that  the  probe 
shaft  was  interferring  with  the  flow  over  the  sensors  and 
might  in  fact  be  generating  a region  of  separated  flow,  it 
was  theorized  that  the  directional  calibration  functions  f(0,  <}>) 
might  not  be  independent  of  velocity.  To  verify  this  suppo- 


velocity  (to  more  closely  match  the  velocities  in  the  ex- 
periment) and  new  directional  calibration  data  were  obtained. 
The  angles  used  in  this  second  test  are  not  the  same  as  in 
the  first  test  because  preliminary  calculations  using  the 
first  calibration  indicated  that  <j>  was  almost  always  negative 
in  the  higher  velocity  regions  of  the  flow  downstream  of  the 
rotor.  The  values  obtained  for  the  function  f(9,  <t>)  at  this 
velocity  (159  ft/sec)  are  listed  in  Table  IV.  The  same 
trends  which  were  observed  at  the  lower  velocity  are  still 
apparent  but  the  values  are  somewhat  higher  overall  indi- 
cating that  the  interference  of  the  probe  shaft  has  dimin- 
ished. A linear  interpolation  was  applied  so  that 

V - vL 

f (0,  4>)  = £l(8,  40  + VH  - VL  - [fH(0,  40  - 

f L ( e , 40]  (22) 

where  f,  and  fu  are  the  values  of  the  function  as  determined 
from  the  low  velocity  and  the  high  velocity  calibrations, 
respectively,  and  and  are  the  velocities  at  which 
those  calibrations  were  done.  The  velocity  V was  deter- 
mined from  the  previous  iteration.  Using  this  procedure  in 
the  test  case  resulted  in  reasonable  agreement  between  the 
calculated  axial  and  circumferential  velocities  and  the 

actual  values.  A strong  radially  inward  flow  was  still 

« 

indicated  and  this  was  known  not  to  exist.  It  was  observed 
that  increasing  the  value  of  f(0,  4>)  for  the  first  sensor 
would  tend  to  bring  the  radial  velocity  into  better  agree- 
ment; consequently,  this  function  was  modified  to  be 

7i(0,  40  - fife,  4>  ) + Af 


where  fi  is  the  value  calculated  from  Eq . 22,  Af  is  the  in- 
crement added,  and  Fi  is  the  value  used  in  Eq . 14.  The 
calculation  of  velocity  components  was  carried  out  with 


■•••  • • < i * • *♦  • •#*  - 
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various  Af’s  until  one  was  found  which  gave  zero  radial  flow. 

At  this  point,  the  radial  and  circumferential  velocity  components 
were  in  close  agreement  with  the  expected  values  of  Vj  and  v2 
as  is  indicated  by  the  listing  in  Table  V.  The  axial 
velocity,  V3,  does  not  match  but,  as  was  previously  mentioned, 
the  pneumatic  instrumentation  used  to  make  this  particular 
measurement  is  not  very  accurate.  The  source  of  the  error 
in  f is  presently  unknown  and  if  there  had  not  been  a known 
flow  condition  included  in  the  experimental  data,  it  would 
have  been  impossible  to  evaluate  it. 

Results 

The  velocity  components  for  the  run  in  which  the  rotor 
blades  were  removed  are  now  a good  match  to  the  expected 
values.  However,  the  radial  component  was  forced  to  be 
zero  and  the  method  used  to  accomplish  this  also  affected 
the  other  two  components.  In  order  to  check  the  validity 
of  the  method,  the  calibration  constants  were  used  to  cal- 
culate the  velocity  components  at  various  gapwise  locations 
for  an  undistorted  flow  condition  resulting  in  light  blade 
loading.  The  velocity  components  V] , v2 , and  v3  are  plotted 
in  Figs.  9 through  11,  respectively,  as  a function  of  their 
gapwise  position.  These  velocity  components  may  be  used  to 
calculate  the  relative  velocity  downstream  of  the  rotor, 

W2 , and  the  exit  flow  angle,  g2 . 

W2  = Vi2  + v22  + v32  (24) 

02  = tan  ( V 3 / v 2 ) (25) 

A radial  flow  angle,  ip , may  be  defined  as 

ip  = tan  (Vl/v3)  (26) 

The  relative  velocity,  W2 , is  normalized  by  its  gap  averaged 
value,  W2  , and  plotted  in  Fig.  12.  The  flow  angles  02  and  ip 
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are  plotted  in  Figs.  13  and  14,  respectively.  These  para- 
meters were  integrated  across  the  gap  to  obtain  their  gap- 
averaged  values  which  are  presented  in  Table  VI.  The  velocity 
components  and  flow  angles  obtained  from  pneumatic  measure- 
ments in  the  stationary  frame  are  also  available  in  this 
table  for  comparison.  The  axial  velocity  shows  extremely 
good  agreement  as  does  the  radial  velocity  and  radial  flow 
angle.  The  exit  flow  angle,  62 > does  not  show  this  good 
agreement.  Referring  to  Fig.  13,  one  observes  that  large 
fluctuations  in  exit  flow  angle  are  indicated.  Admittedly, 
the  probe  used  in  the  present  experiment  should  not  be  ex- 
pected to  yield  accurate  measurements  in  regions  which  have 
large  velocity  gradients  (such  as  the  blade  wake).  However, 
even  if  the  magnitude  of  these  fluctuations  is  somewhat  in 
error,  there  is  still  ample  evidence  that  the  flow  angle 
fluctuations  are  of  sufficient  intensity  to  cause  the  pneumatic 
instrumentation  to  give  erroneous  flow  angle  indications. 

Upon  this  evidence,  the  use  of  pneumatic  instrumentation  to 
measure  flow  angles  near  the  rotor  exit  (both  the  rotating 
triaxial  probe  and  the  stationary  pneumatic  wedge  probe  were 
located  approximately  .05  chord  lengths  downstream  of  the 
rotor  exit  plane)  is  highly  questionable.  The  overall 
trends  tend  to  indicate  that  the  correction  made  to  the 
directional  calibration  is  valid  because  good  agreement 
between  the  triaxial  probe  measurements  and  those  made  with 
pneumatic  instrumentation  was  obtained  for  the  empty  tunnel 
(rotor  blades  removed)  and  for  the  lightly  loaded  flow 
condition  except  for  the  variations  in  the  exit  flow  angle,  S2  . 

Conclusions 

Virtually  all  of  the  difficulties  encountered  during 
the  calibration  process  and  in  applying  the  calibration  data 
to  convert  anemometer  output  voltages  to  velocities  may  be 
attributed  to  the  probe  geometry.  It  appears  that  one 
should  avoid  having  the  probe  shaft  located  behind  the  sensors 





at  all  costs.  This  limitation  may  render  it  impossible  to 
mount  a triaxial  velocity  sensor  and  a total  pressure  trans- 
ducer on  the  same  shaft.  Furthermore,  the  relatively  wide 
spacing  of  the  three  sensors  on  this  particular  probe  leads 
to  erroneous  flow  angle  measurements  in  regions  having  strong 
velocity  gradients.  One  of  the  areas  of  prime  interest,  the 
flow  in  the  wake  of  a rotor  blade,  has  these  strong  velocity 
gradients.  Overall,  one  must  conclude  that  the  geometry  of 
the  probe  used  in  the  present  experiment  is  rather  poor. 

The  constant  temperature  anemometer  circuits  which 
were  designed  in-house  especially  for  this  experiment  per- 
formed flawlessly.  There  was  initially  some  concern  that 
the  centrifugal  loads  to  which  the  units  would  be  subjected 
in  this  application  might  lead  to  mechanical  failure  as 
had  been  the  case  with  some  commercial  units  which  had  pre- 
viously been  installed.  These  fears  proved  to  be  unfounded 
and  the  units  did  not  even  require  adjustment  after  the 
initial  installation.  Therefore,  it  is  believed  that  this 
circuit  is  well  suited  to  applications  which  require  the 
tolerence  of  high  acceleration  loads  and  do  not  permit  easy 
adjustment  of  the  anemometer  once  installed. 

In  spite  of  the  problems  introduced  by  the  geometry  of 
the  probe,  calibration  procedures  have  been  developed  which 
accurately  measure  the  directional  characteristics  of  the 
sensors.  For  this  particular  probe,  these  directional 
characteristics  are  velocity  dependent  and  for  one  sensor 
they  appear  to  involve  some  dependence  which  can  not  be  iden- 
tified and  reproduced  in  the  calibration  tunnel.  The  use  of 
a test  case  in  which  the  velocity  components  were  known  per- 
mits the  directional  characteristics  of  this  sensor  as 
measured  in  the  calibration  tunnel  to  be  corrected  so  that 
velocities  and  flow  angles  may  be  computed  from  the  experi- 
mental data. 
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Recommendations 

The  probe  used  in  this  experiment  required  entirely 
too  much  effort  to  achieve  meaningful  velocity  measure- 
ments. In  particular,  the  directional  characteristics  of  the 
sensors  were  velocity  dependent  and  one  sensor  gave  evidence 
of  being  in  a region  of  separated  flow.  The  measuring 
volume,  the  space  within  which  the  three  velocity  sensors 
are  located,  is  too  large  to  permit  accurate  measurements 
in  flows  having  high  velocity  gradients.  It  is  considered 
imperative  that  a probe  geometry  be  found  which  does  not 
entail  these  difficulties  before  extensive  rotating  frame 
measurements  are  attempted. 

From  the  flow  angle  measurements  made  with  the  triaxial 
probe,  it  is  apparent  that  there  are  severe  flow  angle 
fluctuations  in  the  near  wake  of  the  rotor.  Therefore,  it 
is  certain  that  flow  angle  measurements  made  near  the  rotor 
exit  plane  with  pneumatic  instrumentation  (such  as  wedge 
probes)  will  be  in  error  to  some  extent.  It  is  recommended 
that  the  flow  angle  measurements  in  any  future  experiments 
be  made  with  hot-wire  or  hot-film  probes  so  that  the  instan- 
taneous flow  angles  may  be  accurately  determined  and  a 
proper  average  angle  calculated. 
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Appendix  A 

The  following  pages  are  a listing  of  a program  which 
converts  the  data  in  a table  to  a polynomial  curve  fit  in 
two  varibles.  A least  squares  technique  is  employed  to 
achieve  a certain  degree  of  smoothing  of  the  data.  The 
subroutine  referred  to  as  SIMQ  is  a system  library  routine 
for  solving  simultaneous  equations  of  the  form 

A x = B. 

The  solution  vector  is  returned  in  B and  the  coefficient 
vector  A is  destroyed.  Any  system  routine  which  performs 
the  same  function  could  be  substituted.  In  the  event 
that  one  is  not  available,  a listing  of  SIMQ  is  included. 
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The  polynomial  coefficients  may  be  used  to  reconstruct 
the  function,  f,  at  any  arbitrary  value  of  theta  and  phi  by 
the  algorithm  listed  below.  Since  this  algorithm  will  return 
a result  for  values  of  the  angles  far  outside  the  range  of 
the  original  data  (and  this  result  may  well  be  in  error 
because  the  polynomial  fit  is  not  a good  means  of  extrapo- 
lation) , one  must  be  careful  in  applying  it. 


FUNCTION  FPT(PHI, THETA, N) 
COMMON/ DCOF/A (3 ,5,5) 

FPT=0 . 0 
DO  200  1=1,5 
C=A(N, 5,6-1) 

DO  100  J=2 , 5 

100  C=C*THETA+A(N , 6 - J , 6 - I ) 

200  FPT=FPT*PHI+C 
RETURN 


Table  I.  Response  of  the  Three  Velocity  Sensors  to  Normal  Flow. 
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Values  Obtained  for  the  Function  f(0,  cj))  at  a velocity  of  116  ft/sec. 
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Table  V.  Results  from  Test  with  Rotor  Blades  Removed. 


Parameter 

Pneumatic 

Measurement 

Triax  Probe 
Measurement 

V1 

0 

0.13 

V2 

-7.27 

NO 

o 

00 

1 

V3 

68.70 

76.69 

Table  VI. 

Results  from 

Condition. 

Lightly  Loaded 

Undistorted  Flow 

Parameter 

Pneumatic 

Measurement 

Triax  Probe 
Measurement 

W2 

14  3.5 

155.9 

02 

64.2 

62 . 3 

0.0 

0.5 

c 

129.0 

131.8 

Note:  Velocities  are  in  ft/sec  and  angles  in  degrees. 
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Figure  2.  Rotor  and  Sensor-based  Axis  Systems. 


NOTES:  1.  Op  Amps  are  LM301A  with  30pf  compensation  capacitor. 

2.  Resistors  in  ohms.  1/4  watt  5-6  unless  noted. 

3.  Capacitors  in  microfarads. 

4.  Probe  operating  resistance  is  3 times  control 
resistor . 


Figure  3.  Schematic  Diagram  of  Constant  Temperature 
Anemometer  Circuit  Used  in  LSRR. 
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Figure  6. 
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